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We describe a coupled-chain construction for chiral spin liquids in two-dimensional spin systems. 
Starting from a one-dimensional zigzag spin chain and imposing SU(2) symmetry in the framework 
of non-Abelian bosonization, we first show that our approach faithfully describes the low-energy 
physics of an exactly solvable model with a three-spin interaction. Generalizing the construction 
to the two-dimensional case, we obtain a theory that incorporates the universal properties of the 
chiral spin liquid predicted by Kalmeyer and Laughlin: charge-neutral edge states, gapped spin-1/2 
bulk excitations, and ground state degeneracy on the torus signalling the topological order of this 
quantum state. In addition, we show that the chiral spin liquid phase is more easily stabilized in 
frustrated lattices containing corner-sharing triangles, such as the extended kagome lattice, than in 
the triangular lattice. Our field theoretical approach invites generalizations to more exotic chiral 
spin liquids and may be used to assess the existence of the chiral spin liquid as the ground state of 
specific lattice systems. 
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I. INTRODUCTION 

Understanding the ground states of frustrated quan¬ 
tum spin systems—in which the local energetic con¬ 
straints cannot all be simultaneously satisfied—is a fas¬ 
cinating topic in condensed matter physics^. One of 
the central proposed ground states is Anderson’s res¬ 
onating valence bond state^, a collective spin singlet not 
breaking any symmetry and possessing neutral spin-1/2 
excitations. This idea opened the way for topological 
phases with fractionalized excitations emerging in frus¬ 
trated spin systems^^^^. In 1987, Kalmeyer and Laugh- 
lin^^ proposed a different spin singlet state in the tri¬ 
angular Heisenberg antiferromagnet that breaks time re¬ 
versal and parity symmetries, called the chiral spin liq¬ 
uid (GSL). In 1989, Wen, Zee and Wilczek^^, and also 
Baskaran^^, proposed to use the expectation value of 
the “spin chirality operator” Si • (Sj x Sfc), where i,j,k 
belong to an elementary triangle, as an order parame¬ 
ter for CSLs. Despite preserving spin SU(2) symmetry, 
the Kalmeyer-Laughlin CSL shares basic properties of 
quantum Hall states, such as a bulk gap and chiral edge 
states^^’^®’^®. 

While it was shown later that the CSL is not realized in 
the Heisenberg antiferromagnet on the triangular lattice, 
a few models have been proposed^^“^° for which the CSL 
state is an exact ground state. However, the question 
remained as to whether the CSL can be realized in more 
realistic spin models. Recently, along with related im¬ 
plementations using ultracold atoms in optical lattices^^, 
Bauer et studied a model of a Mott insulator on 
the kagome lattice using exact diagonalization and den¬ 
sity matrix renormalization group (DMRG) and found 
unambiguous evidence for realization of the Kalmeyer- 
Laughlin CSL. The model explicitly includes the three- 
spin interaction • (Sj x Sj,), which is generated by 
an applied magnetic field. Furthermore, He, Sheng and 
Chen^^, as well as Gong, Zhu, and Sheng^^, reported a 


numerical observation of a CSL in an extended spin-1/2 
kagomme Heisenberg model including up to next-next- 
nearest-neighbor interactions. The observation was again 
based on DMRG simulations on cylinder geometries, but 
in this case the spin chirality order emerged from sponta¬ 
neous breaking of time reversal symmetry. Remarkably, 
when second and third-neighbor couplings are small, in¬ 
stead of a CSL one finds^^’^^ a gapped Z 2 spin liquid 
which had been identified in previous studies^®. Quite 
recently, variational Monte Carlo results have confirmed 
that the CSL state is energetically favored in a large re¬ 
gion of the phase diagram of the extended kagome lattice 
and has significant overlap with the exact ground state 
obtained by exact diagonalization^^’^®. 

In order to determine whether the ground state of 
a specific lattice model is a CSL, nonperturbative ap¬ 
proaches not restricted to finite systems are desir¬ 
able. Here we present a new field-theoretic approach 
which captures all universal properties of the Kalmeyer- 
Laughlin state, including fractional quasiparticle exci¬ 
tations and degeneracy on the torus. Our approach 
is based on the “sliding Luttinger liquid” or “coupled- 
wire approach” to the fractional quantum Hall effect 
(FQHE)^®’^°. Similar constructions based on arrays of 
one-dimensional (ID) subsystems have proven powerful 
in the description of exotic quantum Hall states and 
non-Abelian anyons^^^^®, fractional topological insula- 
tors®^^^®, liquids of interacting anyons^^’^^ and purely 
ID systems^®-^®. 

We construct a two-dimensional (2D) CSL from an ar¬ 
ray of antiferromagnetic Heisenberg spin chains. Leaving 
the detailed derivation to the bulk of the paper, here we 
describe the construction pictorially. In the limit where 
the spin chains are decoupled, each chain has gapless spin 
wave excitations moving either to the left {L) or to the 
right (i?). A topologically trivial gapped phase of the 2D 
spin system arises if an energy gap is produced due to 
coupling of L and R movers within the same chain. On 
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the other hand, a topologically nontrivial phase arises 
if the energy gap stems from coupling of the L and R 
modes of neighbohring chains. As we demonstrate later, 
this picture implies the emergence of edge states for a ge¬ 
ometry with open boundaries, consisting of the unpaired 
L and R modes in the spatially separated edge chains. 
Since these edge states are charge neutral, the Hall con¬ 
ductivity vanishes; yet, they are able to conduct heat as 
well as spin currents. 

The bulk Hamiltonian in the topological phase locks 
the L and R spin modes on neighbohring chains into an 
SU(2) symmetric spin singlet state. The theory predicts 
that the elementary excitations carry spin 1/2 and are 
charge neutral; these are the quasiparticles of the CSL. 
Since the excited states in the lattice with an even num¬ 
ber of sites must have integer spin, the spin 1/2 elemen¬ 
tary excitations are fractional and the ground state has 
topological order compatible with filling factor v = \/2 
in the FQHE description of the CSL^^. As compared 
to the electronic FQHE, the extra spin SU(2) symmetry 
implies that quasiparticles and quasiholes—equivalent to 
spin-up and spin-down states—are degenerate. 

The topological nature of the phase, which accounts for 
its long-range entanglement^®, can be tested by placing 
the 2D surface on the torus and counting the ground state 
degeneracy^^. This degeneracy emerges when the oper¬ 
ators that transport an elementary quasiparticle along 
the two non-contractible directions of the torus do not 
commute. In our construction these operators have a 
natural bosonized expression which shows their noncom¬ 
mutativity and the resulting doubly degenerate ground 
state, again consistent with the defining properties of the 
v = ll2 FQHE. 

The most crucial condition for the applicability of our 
approach is that the coupling between L and R spin 
modes of neighbohring chains opens an energy gap. To 
establish the feasibility of this condition, we start by 
analyzing a model of a zigzag chain containing chiral 
three-spin interactions, for which exact results by Frahm 
and Rodenbeck^® provide direct support to our approach. 
This agreement invites the extension to frustrated 2D lat¬ 
tice models, e.g. variants of the triangular and kagome 
lattices, in which recent numerical calculations observed 
signatures of the CSL. The field theory construction sim¬ 
ilarly opens the way for generalizations to more exotic 
chiral spin liquid phases beyond the Kalmeyer-Laughlin 
state. We shall discuss this in the outlook section and 
leave a detailed study for future work. 

The paper is organized as follows. We start in Sec. H 
with one spatial dimension. We first review the ex¬ 
actly solvable lattice model introduced by Frahm and 
Rodenbeck^® and then apply non-Abelian bosonization 
techniques to recover its low-energy physics, forming the 
basis of our wire construction in the simplest context of 
two chains. We show that the spin chirality operator 
opens only a partial gap in the spectrum, leaving out 
two gapless modes which are the seed of the chiral edge 
modes in the 2D case. The 2D construction is done in 


Sec. HI, where a renormalization group analysis is carried 
out to study the competition between the CSL and other 
conventional instabilities. In Sec. IV the properties of 
the ground state obtained in this chain construction are 
discussed, starting with the edge states. Then the quasi¬ 
particles and their creation operator are constructed in 
Sec. IV A, and the algebra leading to the ground state 
degeneracy is described in Sec. IV B. Finally, in Sec. V 
we conclude and discuss future directions. 


II. ONE-DIMENSIONAL CHIRAL SPIN LIQUID 

In this section we use field theory methods to analyze 
a spin-1/2 zigzag model which (i) captures the physics 
of the CSL in one spatial dimension and (ii) is exactly 
solvable"^®. This model corroborates that our field theory 
construction of CSLs, which will become more abstract 
in the next section on 2D generalizations, can indeed de¬ 
scribe concrete lattice realizations. 


A. Spin model 

We analyze a spin-1/2 zigzag chain as shown in Fig. 1, 
described by the Hamiltonian H = Hj + H^. Here 

iL,r = ^[J'S,-S,+i-f JS,-S,+2]. (1) 

For dominating nearest-neighbohr antiferromagnetic ex¬ 
change J' > 0, this system behaves as a single chain 
perturbed by next-nearest-neighbohr coupling J. The 
ID Heisenberg model with J = 0 is exactly solvable by 
Bethe ansatz and the ground state is in a critical phase 
with quasi-long-range order^®. It is known that upon 
inclusion of a small next-nearest-neighbohr coupling J 
the system remains critical, till an energy gap opens for 
J! J' > 0.241167®®“®^. On the other hand, in the limit 
J J' ovl which we will focus, the system can be thought 
of as two chains weakly coupled by the zigzag term J'. To 
force the system into a chiral spin state, we add terms 
breaking parity and time reversal symmetry explicitly 
(but preserving the SU(2) symmetry) 

Hx = fE[Sv(S2,+ixS2,_i) 

j 

-I-S 2 J + I • (S2j X S2J+2)]- (2) 

In both terms in Eq. (2) the spin operators appear clock¬ 
wise in the triple product (as read from left to right) with 
respect to every elementary triangle in Fig. 1. Thus, 
this interaction favours uniform spin chirality. Note that 
the system is not invariant under translation by one site 
j ^ j + 1, but only under j ^ j + 2. 

It is worth mentioning that such chiral three-spin in¬ 
teractions arise naturally in the Hubbard model in the 
presence of a magnetic flux. Following Ref. [22] , consider 
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J 2j-l 2j + 1 

FIG. 1: Spin-1/2 zigzag chain with nearest-neighbor exchange 
coupling J', next-nearest-neighbor coupling J and three-spin 
interaction x- In each triangle the spins are coupled via the 
spin chirality operator, with the order in the triple product 
as indicated by the arrows (see Eq. (2)). 



spin-1/2 electrons hopping on the same zigzag lattice 
with nearest-neighbor hopping amplitude ti and next- 
nearest-neighbohr hopping t 2 - At half filling and for 
strong on-site Hubbbard interaction U we ob¬ 

tain the usual exchange coupling in Eq. (1) with J' = 
J = At^/U. Adding a magnetic flux 0 < $ < tt 
through each triangle breaks time reversal symmetry and 
gives rise to an interaction involving the spin chirality op- 
erator as in Eq. (2), with x ~ sin(<i)). Higher orders 

in t/U, which are required to describe weak Mott insula¬ 
tors, tends to enhance the ratio 

B. Exact spin-wave spectrum 

The Hamiltonian H = Hj + is integrable if one 
parametrizes the three coupling constants as^® 

J' = 2(1 — k), j = k, X = 2-\/ /c(l — k). (3) 


FIG. 2: Spinon dispersion for the exactly solvable model of 
Eqs. (1) and (2) with couplings constrained as in Eq. (3). The 
three curves correspond to different values of k; from bottom 
to top: K = 0.99 (solid line), k = 0.97 (dashed line), and 
K — 0.95 (dotted line). 

Hubbard model, where the operator destroys an elec¬ 
tron with spin a on site j. At long distances compared 
to the lattice spacing a = 1, we expand the fermion field 
around the left and right Fermi points k « ±7r/2 and 
introduce chiral fermions 'ipL,R,cr 

c,,, ^ ~ (5) 

The chiral fermions can be subsequently bosonized as 

a = L,R = +,-, (6) 

where ipa,a are chiral bosons that obey the commutation 
relations 

(^) 5 (^ )] — x). (7) 


Varying the parameter k from 0 to 1 interpolates be¬ 
tween a single Heisenberg chain and a pair of decoupled 
chains. The excitation spectrum e{k) of elementary exci¬ 
tations — called spinons — has a closed-form expression 
extracted from the Bethe ansatz solution"^®, and is plot¬ 
ted in Fig. 2 for different values of k near 1. For k = 1 the 
system reduces to two decoupled chains and one observes 
two branches of excitations containing left- (L) and right- 
(i?) moving gapless modes at k mod 27r = 0 and k = ir 
(in units where the lattice spacing a = 1). For arbitrar¬ 
ily small deviation of k from unity, the pair of R and L 
movers at fc = tt acquires an energy gap^® 

A = e{k = tt) = 27rJe~x, (4) 

while the pair at k mod 27r = 0 remains gapless. 

Below we will obtain this behavior using field theory 
methods. It will also be possible to explain the scaling of 
the energy gap with the interchain coupling. We begin by 
setting the notation, starting from the Hubbard model. 


C. Bosonization notation 

We follow the notation of Ref. [53] and for complete¬ 
ness include the main formulas here. We start from the 


We then introduce charge and spin degrees of freedom 
‘Pa,tix) + ‘Pa,iix) 


^a,c{x) = 
Va,s{x) = 


- ‘Pa,i{x) 


( 8 ) 


At half filling, an arbitrarily small U > 0 gaps out the 
charge mode; this happens through the umklapp opera¬ 
tor, whereby two electrons of opposite spin scatter from 
the right to the left Fermi point and vice versa. The 
low-energy properties are then described by the spin dy¬ 
namics. From now on we will omit the spin index s from 
the spin boson, ipa^s —t ‘Pa- 

The expansion of the spin operator = cj ^ 
reads®^ 


S, ^ S(x) ~ Jfl(x) + 3 l { x ) + (-l)"n(x). (9) 


The spin field contains two parts. The uniform part is 
given by the chiral currents J_r,l(x) of the SU(2)i Wess- 
Zumino-Witten (WZW) model with central charge c = 1. 
In Abelian bosonization notation. 


Ja{x) 
Ja{x) = 


V47r 

2tt 


(10) 
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The staggered part of the spin operator can be written chain and link parity for the odd chain) takes x —)■ —x 

and 


n(a;) oc tr[g(a;)cr], (11) 

where cr is the vector of Pauli matrices and 

g(x,T) = g^iz) 0 gjj(z), (12) 

is the matrix field of the WZW model, with components 
gaa'{x,T) = gL,ai^)9R,a'i^)- Here z = VsT + ix and z = 
VgT — ix are complex coordinates in Euclidean spacetime, 
with Vs the velocity of the spin mode. The spinor fields 
dL.R fields have conformal dimensions®^ (j) 0 ) and ( 0 , |), 
respectively; in Abelian bosonization they can be written 

ga{x) = ^ j ■ ( 13 ) 

D. Interchain coupling 


V : 3lr^3i,L, He ^ He, TClo -)■ -Uo, 

£e —> —Eei £o So- (18) 

In the absence of the chiral three-spin interaction 
ix = 0)j the interchain couplings must respect £, V 
and T symmetries, as well as SU(2) invariance. It is 
known®®’®^’®® that in this case the leading perturbations 
to Eq. (15) are all marginal operators. First, even 
the decoupled-chain Hamiltonian is perturbed by the 
marginal “backscattering” ®® operator 

SHbs = 27 rus 7 bs ^ / dxJi^L ■ J/.i?., ( 19 ) 

/—e,o 

with 7 bs < 0 ; in addition 7 bs = 0{1) since it stems from 
the intrachain exchange coupling J. Second, there is the 
interchain current coupling 


We now turn to the coupling between two Heisenberg 
chains in the zigzag geometry and the resulting phases. 
This has been the subject of extensive theoretical work, 
see for example Refs. [52,56,57]. Here we focus on the 
role of the spin chirality operator H-^. 

In the continuum limit for two weakly coupled chains 
(J',X J), we write the spin operator in even (e) and 

odd (o) chains as 


dHg — ‘I'XVsg J dx(Je,i? * T ' *le,L); (^9) 

where g ~ 0{J'/J) is a dimensionless coupling constant. 
Finally, there is the “twist” operator®^’®® 

= 27rus7tw / dxng • dx-O-o, ( 21 ) 


S2i Se(a:) ~ 3e,Lix) + 3e,Rix) + (-l)“ne(a;), 
S 2.+1 ^ So(x) ~ JoAx) + ^oAx) + (-l)"n„(a;), (14) 

and write the free Hamiltonian in Sugawara form 
2ttVs 


Ho = 


l—e,o 


^ J dxi3lj, + 3fA. (15) 


The Hamiltonian in Eq. (15) describes two pairs of gap¬ 
less right- and left-moving bosonic fields, each pair prop¬ 
agating in one chain. 

To analyze the perturbations to Hamiltonian (15), 
let us discuss the operator content of the theory 
and the symmetries of the lattice model. All lo¬ 
cal operators in the WZW model can be expressed 
in terms of the dimension -1 chiral currents 3i^a{x), 
the dimension- 1/2 staggered magnetization n;(a::) and 
the dimension-1/2 (SU(2) scalar) dimerization operator 
£i{x) oc tr[g;(a;)]®®’®®. These operators transform under 
translation x ^ x-\-l [i.e. j —>■ j + 2 in the zigzag chain) 
in the form 


which carries nonzero conformal spin. The dimensionless 
coupling constant is 7tw 0{J'/J). One can then ana¬ 
lyze the renormalization group (RG) flow of the marginal 
coupling constants. Here we have neglected the marginal 
current-current coupling of the form Jo,fi+(R —>■ L), 

which does not renormalize to one-loop order®^. For a 
single Heisenberg chain, the intrachain coupling 75 ^ < 0 
is marginally irrelevant. On the other hand, for antifer¬ 
romagnetic interchain coupling J' > 0, both g and 7 tw 
flow to strong coupling, but g reaches strong coupling 
first®^. In this case the zigzag chain is in a topologically 
trivial dimerized phase in which both pairs of right and 
left movers within each chain are gapped out. 

To recover the spectrum discussed in Sec. HB, we now 
consider the effects of a nonzero three-spin interaction. In 
this case we must allow for perturbations that are odd 
under V and T, but invariant under the product V oT- 
We find that there is only one new marginal perturbation 
to the decoupled-chain Hamiltonian (and still no relevant 
perturbations). The latter can be obtained by taking the 
continuum limit in the spin chirality operator in Eq. (2): 


C- : J;,a Ji,a) h; —>■ —n/, £/ —^ —£/. (16) 

Time reversal T acts as follows: 

7” : O —Jz,L, Hz ~riz, £i —)■ £[. (17) 

Reflection V about an axis perpendicular to the chains 
that goes through an even site (site parity for the even 


Hx ^ ^ J dx {Seix) ■ [Soix + 1) X So{x)] 

+So(a:) • [Se(a:) X Se(:r + 1)]}. (22) 

We now substitute the mode expansion Eq. (14) into Hx- 
Since we have two field operators appearing at nearby 
positions, we must take their operator product expansion 
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(OPE)®*^ on the same chain. We can use the OPE for 
chiral currents^® 


^ab 
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~ g^ 2^2 + ■ • ■ )(23) 

where is the Levi-Civita antisymmetric tensor. This 
leads to 

Si{x + 1) X Si{x) ~ -[J/,l(x) - Ji,i?(x)] + ... (24) 

TT 



e = iii(A„/A) 


as the leading term. Keeping only nonoscillating terms 
in Eq. (22), we obtain the marginal perturbation 

~ (Je,i? ' Jo,L Je,L ' Jo,i?)5 (^b) 


FIG. 3: RG flow of marginal coupling constants A+ (solid 
line), A- (dashed line) and ytw (dot-dashed line) in the zigzag 
chain, according to Eq. (32). Here we set the initial conditions 
X(0) = 0.1, g(0) = 7tw(0) = 0.05 and yts = -0.5. The 
nonuniversal prefactor is set to C = 1. 


with dimensionless coupling constant x ~ ^{xlJ)- As 
expected, the operator in Eq. (25) is odd under both 
V and T. We remark that we have also considered the 
OPE of the n(x) field with itself and with the chiral cur¬ 
rents, but found no additional perturbations at the level 
of marginal operators. 

Using the OPEs®®’®^ among the fields n; and ei , 
we derive a set of coupled RG equations for the marginal 
coupling constants in the presence of the chirality oper¬ 
ator: 


dXhs 

M 

dg 

M 

dx 

d£ 

dytw 


7te> 

(26) 

g^+je+2T,^C^xL, 

(27) 

25X, 

(28) 

{g + 7bs)7tw, 

(29) 


where (7 is a nonuniversal prefactor of order unity ap¬ 
pearing in the OPE of n(x) and d£ = (iln(Ao/A) with A 
the ultraviolet momentum cutoff. 

Let us now discuss the RG flow. First note that 7bs 
starts off with a bare value of order 1, but is marginally 
irrelevant for 7bs <0. To analyze the remaining equa¬ 
tions, let us assume x > 0 without loss of generality. It 
is convenient to define 


A± = X±5- 


(30) 


The combination of the g and x marginal perturbations 
can be written in the form 

SHg+6H^ = 2TrVsX+ J dxJe,R ■ 3 o,l 

+2TrVsX- J dx3e,L ■ ^o,R- (31) 


The RG equations Eqs. (27)-(29) become 

A^ + 2n^C^xL, 

-A^ - 27r^C^xL 

(^^^+ 7 bs) 7 tw. (32) 

We are interested in the regime x > g, as follows from 
Eq. (3) with 1 — k 1. Physically, a sizeable three- 
spin interaction x can be generated from virtual electron 
hoppings in a Mott insulator in the vicinity of the metal- 
insulator transition^^’®^. In this case A+ > 0 flows to 
strong coupling while A_ < 0 flows to zero. Notice that 
the twist operator contributes to enhancing this trend. 

Figure 3 shows a typical example of RG flow for a given 
choice of bare coupling constants. The important point 
is that A+ reaches strong coupling {i.e. becomes of order 
I) first. This behavior is robust for a wide range of initial 
values in the regime X> g- This implies that, in order to 
understand the properties of the low-energy fixed point, 
we can analyze the effects of large A+ while dropping the 
other competing marginal operators. According to Eq. 
(31), the limit of coupling A+ —)■ oo gaps out right movers 
in the even chain and left movers in the odd chain, but 
leaves the pair of modes Je,L, Jo.i?, gapless. Furthermore, 
for small bare X+{i = 0) the gap in the Jo,l pair 
(see Eq. (4)) is exponentially small since the perturbation 
is only marginally relevant. For the opposite chirality, 
X < 0, the same picture holds upon interchanging A+ O 
A_ and reversing the pair of gapless modes. 

The picture we have just described agrees with the 
low-energy spectrum for the integrable model discussed 
in Sec. IIB if we identify the gapped modes with the 
excitations at A: « tt in Fig. 2. Moreover, the field the¬ 
ory analysis shows that the ID chiral spin liquid phase 
is generic {i.e. does not depend on the fine tuning of 
coupling constants in Eq. (3)) and is governed by the A+ 


dX+ 

hf 

dX- 

hf 

d7tw 

d£ 











6 


operator, which gaps out a pair of left- and right-moving 
spin currents in neighboring chains. This result suggests 
a generalization to two dimensions, which we shall dis¬ 
cuss in the next section. 


III. TWO-DIMENSIONAL CHIRAL SPIN 
LIQUID 

In the previous section we showed that the combination 
of the time-reversal-even interchain coupling 5Hg and 
time-reversal-odd 5H^ can lead to a phase with gapless 
chiral modes propagating in different legs of the zigzag 
chain. We now extend the argument to frustrated 2D 
lattices built out of weakly coupled chains, such as the 
spatially anisotropic triangular lattice depicted in Fig. 4. 
The Hamiltonian is of the form H — Hj -|- , where 



J 


FIG. 4: Spatially anisotropic triangular lattice with exchange 
couplings J ^ J' and chiral three-spin interaction x- 

with 7bs ~ Cl(l). The leading perturbations cou¬ 
pling first-neighbor chains are the generalizations of the 
marginal operators discussed in Sec. IID: 


= (33) 

ij 

contains the exchange couplings Jij = J for i,j nearest- 
neighbor sites along horizontal links of the lattice (be. 
within the same chain) and Jij = J' J for i, j nearest- 
neighbor sites along diagonal links (be. in neighboring 
chains). In addition, the Hamiltonian contains three-spin 
operators with uniform chirality x, 

^x = f E S.-(S,xSfc), (34) 

i,j,kGA 

with the indices i,j,k appearing clockwise in each ele¬ 
mentary triangle. 

The phase diagram of the spatially anisotropic trian¬ 
gular lattice has been studied in great detail in the case 
of time-reversal-invariant interchain coupling®®’®^. The 
starting point is a collection of N decoupled antiferro¬ 
magnetic Heisenberg chains, each of which is described 
by an SU(2)i WZW theory. The spin operator at posi¬ 
tion X in chain I is represented by 


Si{x) ~ Ji,l{x) + Ji.r{x) + {-iyni{x), (35) 

where now I = 1,..., N is the chain (or leg) index. The 
free Hamiltonian for decoupled chains reads 



Now consider the perturbations that are allowed by 
symmetry. Besides SU(2), translation C and time rever¬ 
sal T, it is important to take into account the V symme¬ 
try defined in Eq. (18), which for an arbitrary number 
of chains takes x —> —x and 


'■ ■F)' J/,l 5 n; —^ (—l)*ni, ei —>■ (—l)*~'’^£i. (37) 

As before, the intrachain backscattering process gives rise 
to the marginal perturbation 

SHhs = 27r-Cs7bs E dxJi^R - Ji^L, (38) 

i 


SHg = 27rusgE / dx {Ji^r ■ 3i+i^l + RL)l39) 
i '' 

= 27rUs7twE(~^)^ / dxni ■ (40) 

i 

SH^ = 27rUsxE (J/,L • J/+x_R - i? O L)(41) 

i 

where we have included the T-breaking perturbation that 
stems from the three-spin interaction x- The perturba¬ 
tive RG equations for these marginal coupling constants 
are the same as in the ID case, Eqs. (26) through (29). 

However, there is an important difference between the 
ID and 2D cases. As discussed by Starykh and Balents®®, 
for N > 2 there appear two strongly relevant (dimension- 
1 ) perturbations that are allowed by symmetry and cou¬ 
ple next-nearest-neighbohr chains: 

(5ib„ = dxni-ni+2, (42) 

i 

5Hs = VsgeR'Y^ / dxeiei+2, (43) 

where and are dimensionless. Note that these oper¬ 
ators respect the V oj' symmetry with V defined in Eq. 
(37). These are the only allowed relevant perturbations 
even in our case where V and T are separately broken 
by the spin chirality operator. The RG equations for gn 
and read 

f^(l-f).„, (44) 

^ + («) 

where we have included the correction due to the 0{1) 
marginal coupling 7bs. As argued in Ref. [59], the effect 
of 7bs < 0 is to enhance the growth of gn] for gn > 0, 
this favors an instability towards a collinear antiferro¬ 
magnetic phase in which the n/(a;) fields are pinned. By 







7 



e = iii(A„/A) 



FIG. 6: (a) Edge-sharing triangles; (b) corner-sharing trian¬ 
gles. Only in the first case is the relevant coupling gener¬ 
ated at order x^- 


FIG. 5: RG flow for the marginal coupling constant A+ (solid 
line) and relevant coupling constant Qn. For the latter, we 
show two curves corresponding to two different initial condi¬ 
tions: (?n(0) = 0.01 (dashed line) and g„(0) = 0.001 (dot- 
dashed line). The other initial values used in this plot are 
X(0) = 0.1, g = 7tw = 0.05 and 7ba = —0.5. The dashed 
line represents a case in which the relevant coupling starts 
off smaller than A+ but reaches strong coupling first, driving 
the system into a phase with long-range magnetic order. The 
dotted line represents the case in which A+ reaches strong 
coupling first and leads to the GSL phase. 

contrast, for 7bs > 0 (which can happen if one adds a suf¬ 
ficiently large intrachain next-nearest-neighbor exchange 
coupling) a dimerization instability driven by the £i{x) 
fields becomes dominant. 

Naively, one would expect either one of the relevant op¬ 
erators gn or gg to overtake the marginal couplings in Eqs. 
(39-41) and govern the low-energy physics, leading to 
more conventional, long-range-ordered phases. However, 
the fate of the system also depends on the bare values of 
the coupling constants, and cases in which a marginal op¬ 
erator reaches strong coupling before relevant ones have 
been discussed in the literature®^’®^’®'*’. In the following 
we shall assume 7bs < 0 and focus our discussion on gn as 
the most relevant operator that competes with A+. For 
instance, Fig. 5 shows the RG flow for and A+ for two 
different values of the bare g„(£ = 0) <C A+(£ = 0). We 
see that which operator reaches strong coupling first is 
a quantitative question, whose answer is sensitive to the 
precise initial conditions of the RG flow. Nevertheless, in 
the following we shall attempt to make statements about 
the typical qualitative behavior depending on how 5 n( 0 ) 
scales with the interchain couplings J',x <C J. 

Let us then estimate the magnitude of the bare gn(0) 
in Fq. (42). Our original lattice model does not contain 
direct coupling between next-nearest-neighbor chains. 
However, since the operator is allowed by symmetry, we 
expect it to be generated by the RG flow at higher orders 
in J', X- For the triangular lattice without the three-spin 
interaction, Starykh and Balents®® found that > 0 is 
generated during the initial stages of the RG only at or¬ 
der {J'/JY. On the other hand, it has been suggested®^ 
that a ferromagnetic gn(0) < 0 is expected from fluctua¬ 
tions at short length scales at order {J'/JY ■ The latter 


is more consistent with existing numerical results which 
have not observed the collinear antiferromagnetic phase, 
but rather incommensurate spiral order®^’®®. 

In our case, the chiral spin interaction provides another 
source of the relevant coupling between second-neighbor 
chains. Assuming that the initial value ^^(O) is set by 
fluctuations at short length scales, we can argue that in 
the triangular lattice gn is generated at order xY ■ The 
argument is based on a perturbative calculation in a real- 
space picture. We proceed along the lines of Refs. [58, 
66]. Let Hq denote the Hamiltonian of decoupled chains 
and jo) be the corresponding ground state with energy 
Eq. We regard in Eq. (34) as a perturbation to H^. 
We define the projectors P = |0)(0| and Q = 1 — P, 
and write [fllo) = for the projection of an arbitrary 
state jflt). One can then derive an eigenvalue equation 
for jlko) in the form 7Leff|4'o) = F'|^o)j with an effective 
Hamiltonian given by^®’®® 

Hgf{ = Ho +PH^l-RQHY~^RH^, (46) 

where R = {E — Hq)~^ is the resolvent operator. To 
second order in perturbation theory, we can approximate 
R Ki Rq = [Eq — Hq)~^ and iLeff « Ho + ff*® 

effective interaction 

= PHxRoHx- ( 47 ) 

Gonsider now the edge-sharing triangles represented in 
Fig. 6(a). The spins Si in the lower chain and S4 in the 
upper chain both interact with the spins S2 and S3 in the 
middle chain via the chiral three-spin interaction. The 
exchange coupling between Si and S4 can be generated 
at second order in using Eq. (47) and projecting out 
the middle chain. We find 

Riff")^x"G+-Si.S 4 , (48) 

where 

pOO 

G+- = -i dt{T+^{t)T+^{0)) (49) 

Jo 

is the zero-frequency retarded Green’s function for the 
two-spin operator 

T2+3" = - S:^S+). 


(50) 















Note that T^g is equivalent to the z component of the 
spin current flowing between sites 2 and 3 (the choice 
of the component is arbitrary due to SU(2) symmetry). 
Alternatively, via Kramers-Kronig relations^® G'^ can 
be expressed in terms of the dynamical structure fac¬ 
tor for the operator T^~ involving nearest-neighbor sites 
in a single Heisenberg chain. Note that T^~ can 
be viewed as the antisymmetric part (a vector related to 
the spin current operator) of the two-spin tensor operator 
T^j, a,b G {cc, y, z}. To our knowledge, only the dynam¬ 
ical structure factor for the scalar part of this tensor, 
J2a ~ I been calculated by exact meth¬ 

ods®^. In any case, the main point is that by taking the 
staggered parts of the spin operators in Eq. (48) , we ob- 
tain the interaction ^ gn^i ■ n ;+2 (where I and I + 2 
denote the chains that contain Si and S4, respectively) 
with Qn of order 

From the above discussion, we conclude that for the tri¬ 
angular lattice with y > J' we should in general expect 
5 n(0) ~ 0(y/J)^. Since the relevant operator is gen¬ 
erated already at this level of perturbation theory and 
gn{(-) grows exponentially fast with £, regardless of the 
sign of gn(0), it seems rather unlikely that the marginal 
coupling A+ will reach strong coupling first and lead to 
the CSL phase in the triangular lattice, and instead one 
expects a more conventional order®®. 

More propitious conditions for stabilizing the CSL are 
found in lattices where the chiral three-spin interaction 
is confined to corner-sharing triangles. Consider a set of 
spins connected as in Fig. 6 (b). Repeating the pertur¬ 
bative analysis at short-length scales, we verify that pro¬ 
jecting out the spin S 3 in the intermediate corner gives 
rise to a coupling oc Y.a,bi^iS 2 “ SlS^){StSl - S^S§). 
Importantly, the operator — S’JS'f involving the 

spins in the lower chain is odd under link parity but 
even under time reversal. Thus, taking the continuum 
limit cannot produce n/(a;) in the chain containing Si 
and S 2 . Therefore, this O(x^) perturbative calculation 
does not generate the relevant coupling gniii ■ n;+ 2 - At 
the same time, second order in J' generates an opera¬ 
tor cx (Si -I- S 2 ) • (S 4 -I- S 5 ). The operator in the lower 
chain Si -I- S 2 is odd under time reversal but even un¬ 
der link parity; thus, in the continuum limit it does not 
generate ti/ at order {J'Y either. While our argument 
is based on a lattice picture, the same conclusion can be 
reached by integrating out fast modes in the initial steps 
of the RG in the continuum limit (as done e.g. in Ref. 
[59]). Note also that the same arguments can be used to 
rule out ( 7 £( 0 ) at order y^,(J')^. We conclude that the 
corner-sharing triangular geometry has a higher degree 
of frustration in the sense that the relevant coupling be¬ 
tween second-neighbor chains is pushed to higher orders 
in y and J'. In this case we expect ^ 0{xl ■ 

In Fig. 7 we show an example of an anisotropic 2D 
lattice constructed by coupling chains with only corner¬ 
sharing triangles. This lattice differs from the anisotropic 
kagome lattice®® by an additional exchange coupling J 3 
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FIG. 7: Spatially anisotropic 2D lattice with corner-sharing 
triangles, derived from the kagome lattice by adding a third- 
neighbor conpling J 3 along the horizontal lines inside the 
hexagons. 

between third neighbors in the horizontal direction across 
the hexagons. Moreover, our model includes the three- 
spin interaction in each triangle. For J', y <C J 3 < J, the 
starting point for our analysis is that all spins on this lat¬ 
tice belong to a Heisenberg chain. But in this case there 
are two types of chains with different site densities, which 
we call dense chains (with coupling J) and dilute chains 
(with coupling J 3 ). It is worth mentioning that another 
motivation for considering the J 3 coupling is that the 
extended kagome lattice studied in recent DMRG simu- 
lations^®^^® required relatively large (spatially isotropic) 
second- and third-neighbor couplings in order to stabilize 
the CSL phase. 

We can directly apply the perturbative argument 
about corner-sharing triangles to show that between 
second-neighbor dense chains the relevant coupling is 
not generated at order y^, (J')^. Meanwhile, the sites on 
the dilute chains form a triangular sublattice; thus, the 
gn,gs couplings between nearest dilute chains are ruled 
out by symmetry. Therefore, gn and g^ must be fourth 
order in y, J'. In this case, we expect the scenario rep¬ 
resented by the dot-dashed line in Fig. 5, i.e. the initial 
values of gn,<?e are so small that the marginally relevant 
coupling A+ reaches strong coupling first and gaps out 
pairs of R and L modes in neighboring chains. 

In the next section we shall study the properties of a 
2D state dominated by the operator A^. 
showing that this is indeed the Kalmeyer-Laughlin CSL. 

IV. TOPOLOGICAL PROPERTIES 

Building on the results of the previous section, in this 
part we will assume that the perturbation 

N-l „ 

6H+= 2TrVsX+'^ / dx J/,L • (51) 

1=1 

is the leading relevant operator and gaps out pairs of chi¬ 
ral currents in first-neighbor chains. This is 
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similar to the case of the zigzag chain discussed in Sec. 
IID, except that now the modes that remain gapless are 
spatially separated edge states, composed of the right- 
moving spin mode in the I = 1 chain and the left-moving 
spin mode in the I = N chain. This is a concrete realiza¬ 
tion of the idea of merging triangular puddles to form a 
2D topological phase in the network model perspective^^. 
Thus, the low-energy theory of each edge is described by 
a ehiral WZW SU(2)i model. Note that there are still 
symmetry-allowed relevant perturbations that can couple 
these edge modes, but their coupling constants decrease 
exponentially with N and the effect can be neglected in 
the 2D limit. 

Since the edge states do not carry charge, the Hall 
conductivity vanishes. However, at low temperature T 
the chiral edge modes carry an energy current flowing 
counterclockwise around the edge, given by Jq = ffT^, 
where c = 1 is the central charge in this case®®’^°. 

This critical theory of the edge is consistent with the 
properties of the Kalmeyer-Laughlin CSL state. To com¬ 
pletely characterize a 2D topological state, one needs to 
account for the correct bulk physics in addition to the 
edge physics. Below we will discuss the bulk quasiparti¬ 
cles and show that they correspond to spin-1/2 anyons. 
The unambiguous signature of a topological state is its 
degeneracy on the torus. We will demonstrate that the 
state dominated by Eq. (51) is doubly degenerate when 
placed on a torus. This is directly linked with the ex¬ 
change statistics of these spin-1/2 quasiparticles, imply¬ 
ing that they are anyons with statistical phase 9 = 7r/2^^. 


A. Bulk quasiparticle excitations 

Quasiparticle (QP) excitations can be constructed us¬ 
ing a semiclassical picture that follows from the strong 
coupling limit of the interchain coupling Eq. (51). This 
strong coupling picture is easily understood using the 
methods of Kane, Mukhopadhyay, and Lubensky^® for 
the FQHE. In this description, each ID chain consists of 
L and R bosonic fields, g}i,L/R, subject to a cosine per¬ 
turbation of the form cos[A((/?;i)], with some 
constant A that depends on the scaling dimension of lo¬ 
cal operators. In the strong coupling limit, the field dif¬ 
ference — ipi+i^L is localized in one of the minima 
of the cosine potential, and QP excitations are solitonic 
solutions corresponding to jumps between adjacent min¬ 
ima. 

It is instructive to begin our construction of QPs using 
the above method, even though it does not display the 
SU(2) symmetry of the Kalmeyer-Laughlin CSL. This 
symmetry gives rise to extra degeneracies and implies 
that the QPs must transform under an irreducible rep¬ 
resentation of SU(2). In fact, this restricts the pos¬ 
sible fractionalization of quantum numbers in systems 
with SU(2) symmetry. As described below, non-Abelian 
bosonization is the natural language to construct the bulk 
QPs with explicit symmetry properties. 


Let us separate the longitudinal and transverse parts of 
the interchain interaction, J; l ■ Ji-i-i r = i + , 

where 


Of+i = JhJi+i,R. (52) 

+ (53) 

Here 1 + 1/2 represents the link between chains I and 

I + 1. We use Eq. (10) to represent the components of 
spin currents in terms of bosonic fields. The transverse 
part of the interchain coupling yields 

Oz^i = ^ cos[V4//((z3z,l - +1+1,r)]- (54) 

Upon flowing to strong coupling, this operator pins the 
field difference +i,l~+i+i,r to the minimum of the cosine 
potential; hence in the ground state of the gapped phase 

V^{+i,L- +i+i,r) = ‘2Tm + Tr, n G N. (55) 


Consider an excitation in which the argument of the co¬ 
sine in Eq. (54) jumps by ±27r over some hnite region in 
space (the size of the QP, which depends on the energy 
gap). The total spin AS^ accumulated over that region 
is 

AS^ = f + 

i 


\/+K 





(56) 


Thus, QPs have eigenvalues of equal to ±1/2. 

In the above discussion we have ignored the longitudi¬ 
nal operator Of i. We now check that this was legiti- 
mate. Summing over chain index, we can write 

i i 

-{dx+i,R)'^ - {dx+i,L)'^]- (57) 


In the SU(2) symmetric case, the longitudinal and trans¬ 
verse parts of the marginally relevant operator flow to¬ 
gether to strong coupling. As the transverse part Of/i 

1 + 2 

locks the difference tpz.L — +i+i,r in the low-energy limit, 
the first term in Eq. (57) vanishes. The remaining terms 
are equivalent to a renormalization of the spin velocity 
which does not change the qualitative features of the RG 
flow. Thus, the operator Of, i does not affect the strong- 

coupling picture of pinning the fields as in Eq. (55). 

We can also write down the QP creation operator. The 
latter should create a ±27r kink in the field difference 
V^{+i,L — + 1 + 1 ,r) at link 1 + 1/2. Using the commu¬ 
tation relations in Eq. (7), we can easily identify this 
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with the vertex operators (or equivalently 

g^iv^v3i+i,H(x)^ since the chiral bosons are locked together 
in the ground state). Taking the SU(2) symmetry into 
account, we recognize that the vertex operators that cre¬ 
ate QPs with AS'* = -1-1/2 or AS* = ~l/2 are the two 
components of the chiral spinor in the WZW model (c/. 
Eq. (13)): 


(X gi,L{x) 


g^iV^Vl,L{x) 


(58) 


and 6H+ in Eq. (51). Now consider the operator 

(59) 

This is almost the same as above, except that it only in¬ 
volves the total spin in the left-moving chiral sector of the 
Zo-th chain, i.e. Ux = with S;* ^ = J dx J;* ^(a;). 

In abelian bosonization notation, 

Ux = (60) 


The spinor structure of the QP operator makes it explicit 
that it forms a spin-1/2 representation of the SU(2) spin- 
rotational symmetry. This also implies that “particle” 
(AS* = +1/2) and “hole” (AS* = —1/2) excitations 
can be continuously rotated into each other, which of 
course is only possible because the QPs are charge neu¬ 
tral. Furthermore, local physical operators can only cre¬ 
ate pairs of QPs. For instance, the dimerization operator 
ei{x) (X tr[g;(a:)], which is an SU(2) scalar, creates spin 
singlet excitations, whereas the staggered magnetization 
n/(a::) oc tT[crgi{x)] creates triplet excitations. 


B. Topological degeneracy 

The imprint of topological order is a ground state de¬ 
generacy which is sensitive to the topology of the space^^ . 
To establish the ground state degeneracy on the torus, 
it suffices to find two operators Ux , Uy that commute 
with the Hamiltonian but not with each other. This im¬ 
plies that the ground-state manifold must form a repre¬ 
sentation of the algebra obeyed by Ux and Uy which is 
necessarily multidimensional. More generally, the num¬ 
ber of such operators grows with the genus of the sur- 
face^^. In the FQHE case, this algebra is of the form^^ 
UxUy = C/y 1 / 2 , 6 ^’''*/™ with m integer corresponding to 
the filling factor v = 1/m. If we work in the basis of 
the Ux operator and label a ground state by \x), such 
that Ux\x) = x\x), then the state Uy\x) is also a ground 
state but has Ux eigenvalue of ^ x. Only after 

applying Uy m times do we return to the same value of 
X, implying that there are at least m different ground 
states. Recently the coupled-wire approach was used to 
construct the Ux and Uy operators in the FQHE^^. In 
the following we demonstrate this structure for to = 2 in 
our case of a CSL governed by the interaction in Eq. (51). 

One operator that obviously commutes with any lattice 
spin Hamiltonian is , where is the z component 

of the total spin operator in an arbitrary chain I = Iq. 
This operator is either equal to 1 for an even number of 
spins in the chain or to —1 for an odd number. Upon 
coupling with other chains, the total spin of the Zo-th 
chain can only change by an integer; thus, stays 

invariant. As a result, , iJ] = 0. 

Focusing on the low-energy theory, we write our Hamil¬ 
tonian simply as H = Hq -|- 6H+, with Hq in Eq. (36) 


One can check explicitly that [H, Ux] = 0 since H is writ¬ 
ten in terms of and Indeed, using the commu¬ 
tation relations Eq. (7) and the expressions in Eq. (10), 
we verify that ^ changes S';* by ±1, so it does not 
affect Ux- 

We can then work in the basis of Ux- We label states 
by the eigenvalue of S;* ^, which can be split into integer 
and fractional parts, |S;* = \Af + f), where Af G N 

and / = S;* ^ mod 1. Since the eigenvalues of Ux do not 
depend on JV but only on /, we denote these states only 
by the fractional part |/). 

If Ux commutes with all components of J/,o!, how can 
one find a physical operator that does not commute with 
Ux"! The answer is that the theory allows for physical 
operators which do not appear in H but change the spin 
in a chiral sector of a given chain by a fractional value. 
In the WZW model, one such operator is the staggered 
magnetization n/(x) = tT[gi{x)(T]. In abelian bosoniza¬ 
tion, we can write the z component of ri; as 

nf(x) oc sm{y/TT[ipi^L{x) - ipi,Rix)]} 

= h(®)] _|_ (6X) 

This is a dimension-1/2 vertex operator, clearly distinct 
from the dimension-1 chiral currents. Using the commu¬ 
tation relations Eq. (7), we can verify that nf{x) changes 
the eigenvalue of S';* ^ by ±1/2. Thus, the fractional 
part / that labels the eigenstates changes by 1/2, i.e. 
n/(x) switches between the two sectors with |/ = 0) and 

l/ = l/2). 

We can write down a linear combination of nf(x) and 
ei{x) oc cos{^/^^[lpl^L{x) — ipi,ii{x)]} to construct an oper¬ 
ator that commutes with H but does not commute with 
Ux. Consider 




JJ" ^iU^[vi,LUo)-Vi,RUo)] 
I 


(62) 


where xq is an arbitrary position in the chain direction. 
Note that Uy does not change the fractional part of the 
total S;* for any given chain, but it changes the fractional 
part of each chiral sector S/^ separately. The algebra of 
Ux,Uy can be obtained using and 

Eq. (7). Since 

[i^/^T / dxdx(pio,L{x),V^'^ipi,L{xo)] = (63) 

i 
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FIG. 8: Chiral spin liquid on a torus. The operators Ux and 
Uy transport a spin-1/2 quasiparticle along the x direction 
(parallel to the chains) and y direction (perpendicular to the 
chains), respectively. 


we find 


UxUy = -UyUx. (64) 

Rather than demanding [H, Uy] = 0, it is actually suffi¬ 
cient to show that Uy does not change the energy starting 
from any ground state. As long as the system is closed 
into a torus in the y direction and the fields are locked in 
the ground-state manifold according to Eq. (55), we can 
write 


Uy = n ^iy/7T[(pi + i^L{xo)-ipi^L{xQ)] 

I 


(65) 


But in the phase dominated by the relevant perturbation 
Eq. (51) the difference appearing in the exponential in 
Eq. (65) is just a constant; thus, this operator acts as 
a constant in the ground-state manifold. It follows from 
Eq. (64) that 1/ = 1/2) oc 17^1/ = 0) is a ground state 
orthogonal to \f = 0), and our CSL state has the same 
topological degeneracy on the torus as the ^ = 1/2 FQHE 
state. 

The operators U^ and Uy have the physical interpre¬ 
tation of transporting QPs around the torus in the x or 
y directions, respectively(see Fig. 8). Moreover, it 
was demonstrated in Ref. [71] that for anyons with statis¬ 
tical phase 9 the commutation relation between Ux and 
Uy is UxUy = e^^^UyUx- lu our case, this implies that 
the spin-1/2 QPs in the CSL are anyons with statistical 
phase 0 = 7r/2 (ie. semions). 

Having established the correspondence between the 
edge states, quasiparticle properties and topological de¬ 
generacy, we conclude that the state described in this 
section is equivalent to the Kalmeyer-Laughlin CSL. 


V. SUMMARY AND OUTLOOK 

We presented a coupled-chain construction of chiral 
spin liquids. It was first applied as a long-wavelength de¬ 


scription to an exactly solvable lattice model in one di¬ 
mension, and then generalized to two dimensions. In the 
latter case our formulation assumed a dominant relevant 
interchain coupling given by Eq. (51), which stems from a 
chiral three-spin interaction. This formulation yields all 
the universal properties of the Kalmeyer-Laughlin CLS 
state, suggesting their equivalence. 

The energy gap A of the CSL state is exponentially 
small in the chiral interchain coupling x, ie. A oc Je~5, 
reflecting the nonperturbative nature of the present ap¬ 
proach. On the other hand, the smallness of the gap 
raises the concern that there could be other competing 
instabilities that might dominate the low-energy physics. 
By means of a careful renormalization group analysis, we 
found that the interaction responsible for stabilizing the 
CSL can reach strong coupling before other relevant per¬ 
turbations with parametrically small coefficients. This 
scenario is expected in the kagome lattice rather than the 
triangular lattice. While we have considered the spatially 
anisotropic limit of weakly coupled Heisenberg chains, we 
expect the phase discussed here to be adiabatically con¬ 
nected with the CSL observed in recent numerical work 
on the isotropic extended kagome lattice^^“^®77.28^ 

We have not explored here the possibility of sponta¬ 
neous breaking of time reversal symmetry giving rise to 
the spin chirality order parameter. This possibility can 
in principle be investigated using a mean-field decou¬ 
pling of time-reversal-invariant interchain interactions. 
Another interesting question is whether one can use the 
present coupled-chain approach to construct a Z 2 quan¬ 
tum spin liquid^^, which was shown to be stabilized on 
the kagome lattice^®. Many additional possibilities are 
offered by the current approach, including generalizations 
to more exotic chiral spin liquids, for example SU(N) 
CSLs with proposed realizations for ultracold fermionic 
alkaline earth atomsother exotic states can be ob¬ 
tained by starting from SU(2)fe WZW models, which can 
be realized in higher-S' spin chains®^’^^^^^. We leave the 
development of these ideas for future study. 
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